Capillary Filling of Anodized Alumina Nanopore Arrays 
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The filling behavior of a room temperature solvent, perfluoro-methyl-cyclohexane, in ~20 nm 
nanoporous alumina membranes was investigated in-situ with small angle x-ray scattering. Adsorp- 
tion in the pores was controlled reversibly by varying the chemical potential between the sample and 
a liquid reservoir via a thermal offset, AT. The system exhibited a pronounced hysteretic capillary 
filling transition as liquid was condensed into the nanopores. These results are compared with the 
Kelvin-Cohan prediction, with a prediction in which the effect of the van der Waals potential is 
combined with the Derjaguin approximation and also with recent predictions by Cole and Saam. 

PACS numbers: 68.08.Bc, 61.10.Eq, 68.43.h 



Nanoporous materials hold great potential in a wide 
variety of nanotechnology applications such as: DNA 
translocation 0, 0] , nanofluidic transistors Q , templates 
for nanoparticle self-assembly [!, catalysis p, and 
sensors for chemical agents Q; not to mention more 
commonplace applications such as filtering and humid- 
ity sensors ||. The breadth of nanopore research serves 
to underscore the need for a solid physical understanding 
of the evolution of liquids in such systems. Experimen- 
tal studies have been done on porous network systems, 
such as the disordered Vycor [9( and M41S silica ma- 
terials [lo| . Unfortunately, these porous systems have 
very complicated network geometries that make it rather 
difficult to compare the measurements with simple the- 
oretical predictions. Additionally, most of these stud- 
ies are done at low temperatures that are impractical 
for all of the applications described above. Room tem- 
perature measurements on ordered nanoporous systems 
with nearly ideal geometry should provide insight into 
the physical processes governing liquid behavior in more 
general nanoporous systems. 

We describe here in-situ small angle x-ray scat- 
tering (SAXS) experiments of the equilibrium wet- 
ting and capillary condensation of a room temperature 
solvent, perfluoro-methyl-cyclohexane (PFMC), within 
nanoporous alumina (AI2O3, pore diameter ~ 20 nm). 
The anodized alumina system had an ideal geome- 
try described by a parallel arrangement of cylindri- 
cal nanopores with large aspect ratios ~ 1:5,000, see 
Fig. [1] [lll | Experiments were carried out within an en- 
vironmental chamber that allowed precise control of the 
amount of solvent condensed within the pores via changes 
in the chemical potential, A/i, relative to liquid/vapor 




FIG. 1: (a) SEM image shows well ordered hexagonal packing 
of alumina nanopores. (b) SEM of the fractured edge of the 
membrane showing the parallel packing of the nanopore chan- 
nels (seen as vertical lines) normal to the membrane surface. 



Porous alumina membranes were prepared electro- 



coexistence similar to studies of wetting on flat [12j and 
nanostructured surfaces [l^ . Both adsorption and des- 
orption processes were investigated reproducibly via this 
technique. 



chemically using a two-step anodization technique [11 
lH and cleaned as described elsewhere The resul- 
tant alumina membrane consisted of an array of cylin- 
drical parallel pores running the entire thickness of the 
membrane and open on both ends. The macroscopic 
dimensions of the nanoporous membrane were about 
1 cm x 1 cm x 90 microns. Pores formed a 2D local 
hexagonal order (see Fig. Q]) with nearest neighbor dis- 
tances (center to center) of 58±4 nm, and diameter of 
21 ± 5 nm, determined via electron microscopy. After 
cleaning, samples were dried and loaded into a hermeti- 
cally sealed environmental chamber under an atmosphere 
of "ultra-pure" grade N2 . 

The environmental chamber as described elsewhere 0, 
consisted of two concentric cylindrical metal 
chambers that allowed extremely stable (drift of < 
5 mK/hr) temperature control of the sample with a pre- 
cision of ±1 mK. Condensation of liquid solvent into the 
pores was precisely controlled via a positive offset AT 
between the sample temperature, T s and the temper- 
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FIG. 2: SAXS Geometry: (a) Spectrometer setup: A triple- 
bounce monochromator (M) selects Cu K Q i radiation with 
high angular resolution in the horizontal plane. An identical 
triple-bounce analyzer crystal (A) is used after the sample, 
(b) Alumina pores are tilted, ip 3 — 8.5°, in the y—z plane away 
from the parallel condition to take advantage of the enhanced 
horizontal resolution. Scattered intensity is measured as a 
function of 28 in the x — z plane. 



ature, T r , of a liquid PFMC (purum grade, 97% GC) 
reservoir, related to the chemical potential by Ap ss 
H vap AT/T r . [HHEl Here, H v ^ = 33.9 kJ/mol is the 
heat of vaporization of PFMC. Q] For large AT, little 
or no liquid condenses in the pores. As the AT (and thus 
A//) is decreased toward zero, increasing amounts of liq- 
uid condense into the pores until saturation (complete 
volume filling) is obtained. 

During the experiment, the inner and outer cham- 
bers were set at T in = 31.6 °C (T r = 31.9 °C) and 
Tout — 28.5 °C respectively. Although the sample was 
nearly thermally isolated from the inner chamber; for 
high sample temperatures, such as T s = 65.4 °C, the 
large heat loads on the sample increased Tj„ and T r to 
36 °C and 38 °C), respectively. At each AT measured 
however; the stability of T r was better than 50 mK. 

In-situ small angle x-ray scattering (SAXS) data were 
collected using incident monochromatic Cu-K Q i radia- 
tion with high horizontal angular resolution of 0.005° 
(87 /irad) obtained using a triple-bounce channel cut 
Bonse-Hart [l6| Ge(lll) monochromator . High hori- 
zontal detector resolution was obtained by mounting an 
identical triple-bounce Ge(lll) analyzer crystal after the 
sample, in the non-dispersive geometry [lTj (see Fig. [2]). 
Scattered photons were collected with a scintillation de- 
tector. 

The alumina membrane was mounted in a transmis- 
sion geometry with the face of the membrane nearly 
perpendicular to the x-ray beam, aligning the long axis 
of the pores nearly parallel to the beam. In order to 
take advantage of the narrow horizontal angular resolu- 
tion, the membrane was tilted in the vertical direction 



by ip s w 8.5°. This tilt effectively confined the scatter- 
ing to the horizontal plane and minimized wave-guide or 
multiple scattering effects due to grazing angle reflections 
from the pore walls. In this geometry the intensity of the 
(10) powder diffraction peak, due to the 2D hexagonal 
packing arrangement of the co-aligned pores with a well 
defined pore-pore distance, was recorded in-situ as AT 
was varied between 0.1 K up to 27 K for both adsorption 
(cooling) and desorption (heating) paths. 

Immediately prior to collecting the data presented 
here, the pores were cleaned in place by repeated "flush- 
ing" (produced by rapid thermal cycling over the full AT 
range) with PFMC over ten thermal cycles. The inten- 
sity of the (10) diffraction peak gradually increased with 
flushing as compared at the highest AT; then remained 
constant after three cycles. 

As liquid condenses into the pores, replacing N 2 , the 
electron density contrast of the system decreases, eventu- 
ally saturating when the pores are completely filled with 
liquid. This reduction in contrast is measured by a corre- 
sponding decrease in the scattering intensity of the pow- 
der diffraction peaks, I pea k(AT), relative to the dry peak 
intensity, I pea k{dry) (see Fig. |3j). Neglecting absorption, 
the added solvent volume fraction, V(AT), can be deter- 
mined within the small angle approximation 0, [l3| from 
the lowest order peak as follows: 



V(AT) cx ^I P eak{dry) - ^ I peak (AT) 



(1) 



Due to van der Waals (vdW) interactions, the pores are, 
in practice, never "dry" or empty making I pea k(dry) dif- 
ficult to measure. Even at the large AT realized in this 
experiment we still expect monolayer adsorption of mate- 
rial onto the walls of the pores, evidenced by the necessity 
of "flushing" . For small AT, complete volume filling is 
easily observed by the intensity saturation. For the pur- 
poses of this analysis we calculated a theoretical value for 
Ipeak (dry) from the ratio of the form factors squared for 
empty and PFMC filled pores. The form factor, a func- 
tion of the radius, R, liquid film thickness, d, the electron 
densities, ph and ps, for PFMC and AI2O3 respectively, 
and scattering vector q = (2tt/X) sin(26 i ), is as follows: 



F(p L ,Ps,d,R) = 



PlR 



[RJi(qR) -(R- d)Jx(q(R - d))] 



PsR 
Q 



Jl{qR). 



(2) 



The ratio, F(a = 0)/F(a = R) « 0.3, combined with the 
saturated intensity at low AT, provided the necessary 
offset and normalization. 

Values for the fractional volume condensed into the 
pores as a function of AT, derived from the (10) diffrac- 
tion peak intensity for both adsorption and desorption, 
are shown in Fig. [3J X-ray adsorption corrections were 
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AT=(T r /H vap )A M [K] 



FIG. 3: Main: Volume fraction as a function of AT. Adsorp- 
tion, (o), and desorption, (•), data are shown with arrows 
indicating the measurement direction. A vertical line rep- 
resents the Kelvin transition, AT = 1.9 ± 0.2 K. The solid 
curve is a fit of the pre-transition filling within the Derjaguin 
approximation. Horizontal dashed lines indicate the critical 
and metastable volume fractions, V c and V m respectively. In- 
set (LEFT): (1) conformal pore filling with film thickness, d; 
(2) partially filled (metastable condition) pore with menis- 
cus during desorption. Inset (RIGHT): The (10) diffraction 
peak for various AT during the absorption cycle. The arrow 
indicates the intensity decrease with decreasing AT. 



applied to the highest AT points where the heat load 
caused increased vapor pressure in the environmental 
chamber. The system exhibited a pronounced capil- 
lary transition upon filling, occurring at approximately 
AT* = 1.9 ± 0.2 K (at maximum slope). Utilizing the 
simple Kelvin equation for capillary filling: 



AT* = 



An* = 



Tr 



7 

nR 



(3) 



where 7 = 12.5 mN/m is the surface tension of the 
PFMC, and n = 5.1 X 10 3 mol/m 3 is the molar vol- 
ume of the PFMC, this indicates a nanopore radius of 
R = 11.8 ± 1 nm, in agreement with electron microscopy 
measurements. Here vdW interactions have been ig- 
nored. The system also exhibited hysteresis common to 
porous systems [lj| with the desorption transition oc- 
curring at AT* = 2.3 ± 0.2 K. The width of this hys- 
teresis is smaller by a factor of five than the prediction 
by Cohan [20] for a desorption transition at 2AT* due 
to pores emptying by retreating menisci from the pore 
ends. Thus, the simple Cohan argument cannot provide 
a complete description of the hysteresis observed here, 
indicating the necessity of using a more complex theory. 



In their seminal work, Cole and Saam [2l|,|22j put forth 
a theoretical treatment which includes substrate interac- 
tions ignored by the Cohan-Kelvin theory. A brief outline 



of their capillary pore filling theory, including the influ- 
ence of vdW interactions with the pore walls is given 
here. We make an additional assumption, valid for this 
system, that the capillary transition occurs when the 
adsorbed thickness of liquid on the pore walls is thin 
enough, with respect to t he p ore radius, to employ the 
Derjaguin approximation 23j, [2J], allowing a simplified 
analytical treatment. 

The total free energy per unit length of a cylindrically 
symmetric film on the walls of the pore is: 

G(d) = 2Tr~f(R-d)+jiA[iiT{2Rd-d 2 ) + U(d) , (4) 

where U(d) is the vdW potential between the liquid and 
the pore walls. For thin films relative to the pore radius, 
R ^> d, U (d) is approximately that of a liquid film on a 
flat substrate (Derjaguin approximation) [241 ]: 



U(d) 



A 



cir 



12nd 2 



2nR 



(5) 



where A g = 6.3 x 10 J is the effective Hamaker con- 
stant for PFMC on A1 2 3 (2{|. In the limit of vanishing 
film thickness, this term should be replaced by the in- 
terfacial tension term for the pore wall/liquid interface. 
The equilibrium condition is then given by: 



dG{d) 
dN 

A/x(d) 







A c s 7 
6nnd 3 n(R — d) 



(6) 
(7) 



Equation [7] simplifies to the Kelvin equation when vdW 
interactions are ignored. Using Eq. ( [7]) we fit the pre- 
transition filling data for the absorption cycle as a func- 
tion of R yielding a value of R = 11.8 ± 0.2 (see Fig. [3] 
solid curve) in agreement with the Kelvin value above. 
Discrepancies with the fit at large AT are due to the film 
thickness being on the order of the monolayer thickness. 
The capillary transition to filled pores occurs when the 
chemical potential reaches the minimum allowed value, 
since further reduction of the chemical potential can only 
be achieved by complete filling of the pore, d = R, in 
which case the surface tension term vanishes due to the 
vanishing of the interface. This stability or critical con- 
dition for the cylindrically symmetric film is satisfied for 
8Afi(d)/dN > (in the limit d < R): 



A, 



> 



7 



2nd 4 ~ {R-df 



J_ 
R? 



(8) 



This demonstrates that the thin film is stabilized by the 
vdW interactions with the pore wall and that, at the 
transition, d ~ 3.3 nm validating the use of the Derjaguin 
approximation, d <C R. 

For our system with R = 11.8 nm (using the fit value), 
numerical calculations of Cole and Saam predict an on- 
set of instability for volume fraction V c « 0.44 (see Fig. 
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3 in Ref. [22j) with R Q = ^A cS /2j w 1.6 ran for the 
PFMC/AI2O3 system. The estimate for the onset of the 
capillary transition, V c occurs at slightly higher AT than 
expected from the Kelvin and Derjaguin methods (see 
Fig. [3]), but within the uncertainty given above. 

Cole and Saam additionally make a prediction about 
the hysteresis upon desorption. They argue that the hys- 
teresis occurs due to a partially filled pore with menis- 
cus (see Fig. [3j upper insert (2)) being metastable with 
respect to a conformal thin film geometry (see Fig. [3l 
upper insert (1)). The metastable limit is also a func- 
tion of the pore radius, R, and Rq. Cole and Saam's 
numerical analysis predicts a metastable volume fraction 
V m Ri 0.23 for our experimental parameters (see Fig. 3 
in Ref. [13] )• The V m value approximately corresponds 
to the volume fraction where the desorption curve joins 
with the absorption curve, in good agreement with our 
experimental results, as indicated in Fig. [3l It is worth 
noting that this joining occurs at AT ss 2 AT*, where the 
Cohan argument predicts pores to empty by hemispher- 
ical menisci. Unfortunately there is no prediction of the 
AT where the desorption transition initiates, and it not 
possible to translate V c and V m into a hysteresis width. 

The above analysis indicates that for systems with 
ideal capillary geometry, such as the porous anodized alu- 
mina system studied here, the theory of Cole and Saam 
provides a good description of the capillary filling transi- 
tion. For the cylindrical pores described here, the Kelvin 
equation is a reasonable predictor of the filling transition 
as well, though it must become increasingly worse for 
decreasing pore size. The Kelvin-Cohan argument fails 
to qualitatively predict the observed hysteresis for the 
system. Cole and Saam's theory, which includes vdW 
interactions, demonstrates good agreement with the ex- 
perimental results in prediction of the capillary transition 
and does yield information about the observed hystere- 
sis, yet there is not a satisfactory prediction of the ob- 
served width. It is clear from the experimental data that 
the desorption transition occurs at lower AT than either 
prediction. Some of the shift may be due to broadening 
of the transition by pore diameter polydispersity, though 
from Fig. [3] this effect seems insufficient to fully account 
for the shift. Also, large pore aspect ratios rule out the 
influence of end effects and suggests other mechanisms 
for desorption are required to explain this feature of the 
data. 

We are grateful to M. W. Cole for helpful discussions 
and comments. This work was supported by the National 
Science Foundation Grants No. 03-03916 and No. 01- 
24936 



(NIST), Gaithersburg, MD (USA) 
' Current address: School of Chemical and Biological Engi- 
neering, Seoul National University, Seoul (South Korea) 

[1] A. Meller, L. Nivon, and D. Branton. Phys. Rev. Lett. 
86, 3435 (2001) 

[2] Y. Rabin, and M. Tanaka. Phys. Rev. Lett. 94, 148103 
(2005) 

[3] R. Karnik, R. Fan, M. Yue, D. Y. Li, P. D. Yang, and A. 
Majumdar. Nano Lett. 5, 943 (2005) 

[4] K. J. Alvine, O. G. Shpyrko, P. S. Pershan, D. J. Cook- 
son, K. Shin, T. P. Russell, M. Brunnbauer, F. Stellacci, 
and O. Gang. Phys. Rev. B 73, 125412 (2006); K. J. 
Alvine. Thesis, Harard University (2006) 

[5] M. Lahav, T. Sehayek, A. Vaskevich, I. Rubinstein. 
Angew. Chem. Int. Ed. 42, 5576 (2003). 

[6] R. E. Benfield, D. Grandjean, M. Kroll, R. Pugin, T. 
Sawitowski, and G. Schmid. J. Phys. Chem. B 105, 1961 
(2001). 

[7] J. P. Novak, E. S. Snow, E. J. Houser, D. Park, J. L. 
Stepnowski, and R. A. McGill. Appl. Phys. Lett. 83 
4026 (2003). 

[8] O. K. Varghese, D. Gong, M. Paulose, K. G. Ong, C. A. 

Grimes. J. Mater. Res. 17 1162 (2002). 
[9] L. D. Gelb and K. E. Gubbins. Lang. 14, 2097 (1998). 
[10] A. V. Neimark, P. I. Ravikovitch Micropor. and Mesopor. 

Mater. 44, 697 (2001). 
[11] H. Masuda and K. Fukuda. Science 268, 1466 (1995). 
[12] I. M. Tidswell, T. A. Rabedeau, P. S. Pershan, and S. D. 

Kosowsky. Phys. Rev. Lett. 66, 2108 (1991). 
[13] O. Gang, K. J. Alvine, M. Fukuto, P. S. Pershan, C. T. 

Black, and B. M. Ocko. Phys. Rev. Lett. 95, 217801 

(2005) 

[14] K. Shin, H. Xiang, S. I. Moon, T. Kim, T. J. McCarthy, 

and T. P. Russell. Science 306, 76 (2004). 
[15] R. K. Heilmann, M. Fukuto and P. S. Pershan. Phys. 

Rev. B 63, 205405 (2001). 
[16] M. Deutsch. J. Appl. Cryst. 13, 252 (1980) 
[17] J. Als-Nielsen, and D. McMorrow. Elements of Modern 

X-Ray Physics, Wiley, New York, USA, 2001) pg. 200 
[18] W. D. Good, D. R. Douslin, D. W. Scott, A. George, J. 

L. Lacina, J. P. Dawson, G. Waddington. J. Phys. Chem. 

63, 1133 (1959) 
[19] A. W. Adamson, and A. P. Gast. Physical Chemistry of 

Surfaces, 6th ed., Wiley Interscience, New York, USA, 

1997) pg. 665 

[20] L. H. Cohan. J. Am. Chem. Soc. 60, 433 (1938). 
[21] M. W. Cole, and W. F. Saam. Phys. Rev. Lett. 32, 985 
(1974) 

[22] W. F. Saam, and M. W. Cole. Phys. Rev. B 11, 1086 
(1975) 

[23] B. V. Derjaguin, and N. V. Churaev. J. Colloid Interf. 
Sci. 54 157 (1976). 

[24] J. M. Israelachvili. Intermolecular & Surface Forces, 2nd 
ed., (Academic Press, New York, USA, 1992) 

[25] This value for A cft for PFMC on AI2O3 relies on the 
approximation A12 ~ VA11A22 (see Ref. [24]) and the 
values An (PFMC)= 2.8 x 10" 20 J (see Ref. [16]) and 
An (A1 2 3 )= 14 x 10" 20 J (see Ref. [24]). 



* Electronic address: kalvine@post.harvard.edu. Current 
address: National Institute of Standards and Technology 



